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In this paper we consider the quantum phase transition in the Ising model in the presence of 
a transverse field in one, two and three dimensions from a multi-partite entanglement point of 
view. Using exact numerical solutions, we are able to study such systems up to 25 qubits. The 
Meyer- Wallach measure of global entanglement is used to study the critical behavior of this model. 
The transition we consider is between a symmetric GHZ-like state to a paramagnetic product- 
state. We find that global entanglement serves as a good indicator of quantum phase transition 
with interesting scaling behavior. We use finite-size scaling to extract the critical point as well as 
some critical exponents for the one and two dimensional models. Our results indicate that such 
multi-partite measure of global entanglement shows universal features regardless of dimension d. 
Our results also provides evidence that multi-partite entanglement is better suited for the study of 
quantum phase transitions than the much studied bi-partite measures. 
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I. INTRODUCTION 

There has been much work on entanglement in the past 
twenty years Entanglement is a purely quantum 

phenomenon with no classical counterpart. It is thought 
to hold the key to a deeper understanding of the theoret- 
ical aspects of quantum mechanics. From a more prac- 
tical aspect, entanglement is the key ingredient in many 
information processing applications including quantum 
computation and quantum cryptography ,9]. On the 
other hand, there has been much cross-fertilization in 
the fields of condensed-matter physics and quantum in- 
formation theory in recent years |10h18| ] . Here, many tra- 
ditional condensed-matter systems including fermionic 
and bosonic gases, and in particular lattice spin mod- 
els have been investigated in the light of new develop- 
ments in quantum information theory and entanglement 
in particular (Toj . It has been found that entanglement 
plays a crucial role in the low-temperature physics of 
many of these systems, particularly in their ground (zero- 
temperature) state jlftjig ]. A very fruitful avenue along 
these lines has been the relation of entanglement and 
phase transitions in general, and quantum (ground state) 
phase transitions in particular. This is a bit surprising 
since entanglement was originally thought to be some- 
what fragile and thus easily destroyed by fluctuations. 

In a quantum phase transition (QPT) [20], a ther- 
modynamic system described by a Hamiltonian H(X) 
changes its macroscopic phase at the critical value of the 
control parameter A c . In recent studies of many ther- 
modynamic systems exhibiting QPT, in particular quan- 
tum spin 

models dMi], 

it has become clear that the 
onset of transition is accompanied by a marked change 
in the entanglement. Depending on the model, entan- 
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glement could peak, or show discontinuous behavior, or 
show diverging derivatives with scaling behavior at the 
critical point [13, [H|- What is less clear is the exact role 
(or the general mechanism) through which entanglement 
and QPT are related. In such studies, various quantita- 
tively different measures of entanglement have been used. 
Therefore, for example, one would like to know if there 
are universal features in the entanglement of various spin 
model exhibiting QPT? 

Another important feature which is emerging out of 
recent studies of condensed-matter systems from quan- 
tum information perspectives is the need for multi-partite 
measures of entanglement [22hj26| . This, by the way, is 
an example of cross-fertilization referred to earlier. Since 
the root of quantum theory [27j is originally in bi-partite 
systems like Bell states, it has been natural to study 
macroscopic systems using bi-partite measures such as 
the von Neumann entropy or concurrence. In fact, with 
very few exceptions, the general body of the current lit- 
erature has used such bi-partite measures to study many 
particle systems. Although this has been so because of a 
matter of tradition and/or convenience, there is increas- 
ing evidence that such measures are generally inadequate 
to study QPT in condensed- matter systems [28l [20j| . Af- 
ter all, it is natural to use multi-partite entanglement if 
one is to study the role of entanglement in multi-partite 
(many-particle) systems, as important types of entangle- 
ment in such systems (e.g. various n-tangles) may not be 
captured by a bi-partite measure, but would be included 
in an (ideal) multi-partite measure. Additionally, some 
multi-partite measures (as will be discussed in Sec. Ill) 
have thermodynamic properties (e.g. extensivity) which 
make them more suitable for studies of such thermody- 
namic phenomenon as QPT. Another equally important 
shortcoming is that most such studies have been carried 
out for one dimensional (Id) models. Although this is 
perhaps because of computational difficulties, it is cer- 
tainly not well-justified. As is well-known, spatial dimen- 
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sion (d) plays an important role in the physics of ther- 
modynamic systems, phase transitions in particular [30|. 

Here, we propose to study QPT in prototypical trans- 
verse field quantum Ising model using the Meyer and 
Wallach [3l[ measure of global entanglement in one, two, 
as well as three dimensions. Such global entanglement 
measure seems to be well-suited for studies of many par- 
ticle systems [22|, [32j . Since analytic results are usually 
difficult to come up with, numerical results with finite- 
size systems are typically the way to proceed. However, 
solving quantum lattice spin systems numerically is also 
computationally expensive as only a few qubits (spins) 
can be solved exactly and approximation techniques have 
limited success in one dimension and are more limited in 
higher dimensions (33J. Very recently, however, such sys- 
tems have been studied using efficient numerics [34j. 

In this article, we solve the transverse field quantum 
Ising model numerically (exact) for up to 25 qubits in 
one, two and three dimensions. Our main result is that 
global entanglement is a well-suited measure to study 
QPT with some universal features in any dimension. We 
show that global entanglement has interesting scaling 
properties near the critical point. Using finite-size scaling 
arguments, we extract critical points as well as some crit- 
ical exponents for the Id and 2d models consistent with 
previous studies. Due to system-size limits, we are only 
able to study the smallest 3d system and thus cannot 
perform finite-size studies. However, the general shape 
of global entanglement in the 3d model (Fig. 7) indicates 
that our Id and 2d results easily generalize to 3d systems. 
More importantly, our results provide a general frame- 
work for computation of an accessible measure of en- 
tanglement and its relevance to QPT's in many-particle 
thermodynamic systems. 

This paper is structured as follows: in Section II, we 
discuss the multi-dimensional quantum Ising model in 
the presence of a transverse field and its ground state 
properties relevant to our study here. In Section III, we 
discuss some key concepts regarding the Meyer- Wallach 
measure of global entanglement, while our main results 
are presented in Section IV. Our concluding remarks, in- 
cluding suggestions for further work is presented in Sec- 
tion V. 



II. TRANSVERSE FIELD ISING MODEL 

The system under consideration here is the ferromag- 
netic Ising model in a transverse field given by the Hamil- 
tonian: 

N 

<ij> i=l 

where the A = -g where J is the ferromagnetic coupling 
constant, B is the magnetic field, N is the total number of 
spins (qubits) and of and af are the Pauli spin matrices 
in x and z direction at the site i. < ij > means site i and 



j being nearest neighbors on a regular d-dimensional lat- 
tice. We use periodic boundary conditions. This model 
has been extensively studied in Id, but less is known 
about its properties in 2d and 3d. Relevant to our study 
here is the QPT this model exhibits regardless of dimen- 
sion. At zero field this model exhibits ferromagnetic be- 
havior with net magnetization in the x direction while in 
the large field limit, it exhibits a paramagnetic behavior 
where all spins point in the field direction z. The transi- 
tion between these two phases occurs at the critical value 
of A = A c , in the thermodynamic limit. It is well known 
that the ground state is a product state in both these 
limits 35]. In the first limit, the ground state is two- fold 
degenerate, one being a product state of spins pointing 
in the positive x direction |+) = |x; 0)i|x; 0)2---|x; 0)jv 
and the other is |— ) = \x; l)i\x; 1)2. ,.\x; 1) n which is the 
global phase flip of the first one. In the second limit, the 
ground state is a product state of spins pointing in the 
positive z direction |0). Both limits of the ground state 
are product states which are disentangled but there is 
another possibility for the ground state in the first limit 
which arises from linearity of the Schrodinger equation. 
As the |+) and |— ) are solutions for the ground state in 
this limit, the superposition of these degenerate states 
is also another acceptable solution for the ground state 
when the applied field (B) tends to zero. This possi- 
bility is a GHZ-like state which has genuine multi-qubit 
entanglement [36| , 

\GHZ) N = -L(\+) + \-)) (2) 

The possibility of GHZ-like ground state is fascinating 
from a fundamental theoretical point of view as it rep- 
resents a coherent superposition of two macroscopically 
distinct states, and hence, is often called a cat-state. It 
is proven in Ref. [371 ] that in this limit the ground state 
of the Ising model is locally unitary equivalent to an N- 
partite GHZ state, and Ref. [38[ discusses how to prepare 
an Ising chain in a GHZ state using a single global control 
field. It is important to note that such a ground state 
would show zero net magnetization, i.e. < M x >= 0. 
Therefore, such quantity could not be used as an order 
parameter to signal the phase transition under consider- 
ation here. We also note that such ground states have 
recently attracted attention from a symmetry-breaking 
point of view [39j . Here, we study the multi-partite en- 
tanglement properties of such a ground state and its sub- 
sequent transition to a (paramagnetic) product-state as 
a function of A. 



III. GLOBAL ENTANGLEMENT 

Global entanglement, defined by the Meyer- Wallach 
entanglement measure of pure-state [3lj], and henceforth 
denoted by E g i, is a monotone [13], and a very useful 
measure of multi-partite entanglement. As we will show 
briefly, E g i is a measure of total non-local information 
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per particle in a general multi-partite system. There- 
fore, E g i gives an intuitive meaning to multi-partite en- 
tanglement as well as being an experimentally accessible 
measure [UHl. 

Finite amount of information can be attributed to N- 
qubit pure state which is N bit of information according 
to Brukner-Zeilinger operationally invariant information 
measure (43|. This information can be distributed in lo- 
cal as well as non-local form, which is associated with 
entanglement [44j . This information has a complimen- 
tary relation: 

Itotal = Ilocal ~t~ Inon — local • (3) 

The total information is conserved unless transferred 
to environment through decohcrence. The amount of 
information in local form is Ii OC ai = Y]j—-\ h where, 
Ii = 2Trpf — 1 is th e op erationally invariant information 
measure of a qubit [43| , and pi is single particle reduced 
density matrix obtained by tracing over the other par- 
ticles' degrees of freedom. Therefore, according to Eq. 

(3) Inon-iocai = J2iLi 2 (1 ~ Tr Pi) which is distributed 
in different kinds of quantum correlations, the tangles, 
among the system, 

Inon-local = 2 T M2 + — +N n t ...i N i (4) 

ii<i2 ii<— <ijv 

where the first term is referred to as 2-tangle, the next be- 
ing 3-tangle and the last term the N-tangle of the system. 
One can view these tangles as different types of non-local 
information distribution. Therefore, since E g i is the sum 
of single particle linear entropies per unit particle in a 
multi-partite system [3l|, it can be written as: 

e b< = ^ 2 Y. + -+ N Ti ^- ( 5 ) 

il<i2 ii<...<tN 

Therefore, E g i is the average of tangles per particles 
(tvO; without giving detailed knowledge of tangle distri- 
bution among the individual particles. This is much like 
the average energy per particle in an interacting many- 
particle system. The Meyer- Wallach measure was origi- 
nally introduced as a multi-partite entanglement to dis- 
tinguish it from bi-partite entanglement measures like 
entropy of entanglement. But, as shown above, E g i is 
an average quantity and therefore cannot distinguish be- 
tween entangled states which have equal (r) yet differ- 
ent distribution of tangles, like \GHZ) N and \EPR)®% . 
However, E g i can distinguish between GHZ and W states 
since they have different values of (r) , so E g i , like a ther- 
modynamical variable, determines the general amount of 
a property in a quantum system without giving a de- 
tailed knowledge of its sharing among the constituents. 
One expects this property of global entanglement to play 
an important role in studying macroscopic properties of 
multi-partite quantum systems [HI [32j . 

To obtain E g \ in these systems, we have to calculate 
the single particle reduced density matrix, pi. Since we 



use periodic boundary conditions, the reduced density 
matrix is the same for all particles. So E g \ reduced to 
linear entropy of a single particle density matrix, pf. 

E gl = 2(1 - Trpf) (6) 

IV. RESULTS 

Using Eq. (6), we can therefore easily calculate E g i 
exactly for any dimension d, up to the limitations set by 
computational limits of our numerics. 

We start by showing our results for the Id model. Fig- 
ure 1 shows E g i vs A and the inset shows its derivative 
for various system sizes up to N — 24. The general be- 
havior shown here is that of E g i increasing slowly from 
its zero value at A = with a sharp transition to its large 
A value of 1 around A c = 1. The critical point is better 
seen in the derivative (inset) which peaks at the maximal 
value X m (N). As the system size increases, the peak of 
the derivative sharpens and moves closer to the critical 
point A c = 1. 
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FIG. 1: (Color online) Global entanglement as a function of 
A for the Id transverse Ising model. The inset shows the 
derivative and the system sizes used. Increasing N sharpens 
the peak and moves it closer to the critical point. 

The extrapolation to the infinite system size along with 
the convergence to the critical point (inset) is shown in 
Figure 2. As one can see, A m (oo) = 1.01, very close 
to the well-known result of A c = 1, showing that E g i 
is a good indicator of the critical behavior of this model. 
The inset shows that the convergence to the critical point 
is in accordance with |A C — A m | ~ N~ a with exponent 
a = 1.00. 

We next examine the scaling behavior of Egi near the 
critical point. According to scaling ansatz [45(], we have 
dE g i/d\ ~ Q(N»(\ - Am)) where v is the correlation 
length critical exponent, and Q{x) ~ ln(x) is generally 
assumed. As is seen in Figure 3, an acceptable collapse 
occurs for various N's using the scaling ansatz with the 
critical exponent v — 1.06 in line with previous studies 
17], and close to the exact result of v = 1. The inset 
shows the logarithmic divergence of the maximum (peak) 
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FIG. 2: (Color online) Convergence of A m to the critical point 
as TV — !> oo, for the Id transverse Ising model. The y-intercept 
is 1.01. The inset shows the relation |A C — A m | ~ TV -1,00 . 
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FIG. 3: (Color online) Finite-size scaling of global entangle- 
ment for the Id transverse Ising model. The inset shows the 
logarithmic divergence of the value of the derivative at the 
maximal point A m . 



of the derivative of E g i . Hence, the general shape of E g i , 
the logarithmic divergence of its derivative at the critical 
point, along with its consistency with finite-size scaling 
ansatz provides strong evidence for the well-suitedness 
of such measure for the Id Ising model. The question 
now is, if such features also hold for higher dimensional 
models? 

We next turn to the 2d model. Using periodic bound- 
ary conditions wc have been able to study such model for 
up to L 2 = 5 2 = 25 = N qubits. Figures 4, 5 and 6 show 
similar results as that of Figures 1, 2 and 3. We note 
the following: The general shape of the E g i still remains 
(Figure 4) , with a (logarithmic) divergence of the deriva- 
tive at the critical point (insets of Figure 4 and 6). The 
critical point is now identified as A TO (oo) = 0.329(Figure 
5), consistent with recent studies using infinite projected 
entangled-pair state (A c = 0.3268) [34|, as well as quan- 
tum Monte Carlo simulations (A c — 0.3285) [46j. Inter- 
estingly, our simple method obtains a more acceptable 
result than the recent similar multipartite entanglement 
study based on matrix and tensor product states which 



obtained A c = 0.308[24]. We note that the convergence to 
the critical point (inset of Figure 5) is in accordance with 
I A c — A m | ~ L~ a with exponent a — 1.00 being exactly 
the same as the Id case. The difference here is that this 
convergence occurs from above the critical point as op- 
posed to the Id case. The finite-size scaling ansatz is also 
valid (Figure 6), giving the correlation length exponent 
v = 0.51. The inset of Figure 6 shows the logarithmic 
divergence of the derivative at the critical point. 

In Figure 7, we show our result for the 3d version of 
this system for the only system size we are able to study. 
The general behavior of E g i seen in the Id and 2d model 
is clearly seen here for the 3d case as well. While we are 
not able to perform scaling analysis similar to the Id and 
2d models, it seems reasonable to assume that the same 
general behavior carries over to the 3d model. We note 
that X m (L = 2) = 0.26 here, which would understand- 
ably be different from the infinite-size limit, but in the 
right ball-park of A c w 0.2^. 




X 

FIG. 4: (Color online) Global entanglement as a function of 
A for the 2d transverse Ising model. The inset shows the 
derivative and the system sizes used. Increasing system sizes 
sharpens the peak and moves it towards the critical point. 
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FIG. 5: (Color online) Convergence of A m to the critical point 
as L — > oo, for the 2d transverse Ising model. The y-intercept 
is 0.329. The inset shows the relation |A C — A m | ~ L -1 ' 00 . 
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FIG. 6: (Color online) Finite-size scaling of global entangle- 
ment for the 2d transverse Ising model. The inset shows the 
logarithmic divergence of the value of the derivative at the 
maximal point A m . 



FIG. 8: (Color online) Global entanglement (continuous 
curve) and genuine entanglement (dashed curve) for the Id 
(main figure) and 2d (inset) transverse Ising model. The Id 
result is for N = 16 and the 2d result is for a 4 x 4 system. 




FIG. 7: (Color online) Global entanglement and its derivative 
(inset) for the 2x2x2 transverse Ising system, the only 3d 
system we have been able to study. 



Finally, it is worth considering another important form 
of multi-partite entanglement, namely, genuine entangle- 
ment. Genuine entanglement in a many-particle system 
represents the amount of entanglement shared by all par- 
ticles. Therefore, genuine entanglement is equal to the 
./V-tangle, the last term in Eq. (5). One might expect 
that such a term would gradually lose its significance as 
N increases. However, due to the GHZ nature of our 
ground state, this term (iV-tangle) is the dominant term 
in E g i and its dominance increases with increasing A. 
This is shown in Figure 8 for both the Id and 2d models, 
where one can easily see that increasing A, increases the 
share of iV-tangle in E g i . It is also worth noting that the 
structure of genuine entanglement is very similar to an 
order parameter. It is zero on one side of the transition 
and becomes non-zero around the critical point rising to 
its maximum at 1/A = 0. This behavior becomes more 
pronounced as system-sizes get larger, however, we note 
that the transition is not sharp and is in fact "rounded" . 
However, since net magnetization cannot be used as an 
order parameter here, such behavior deserves further at- 



tention. 



V. CONCLUDING REMARKS 

In this paper we have studied the quantum phase tran- 
sition in the transverse field Ising model from a multi- 
partite entanglement point of view on a one, two and 
three dimensional square lattice. Our work is interesting 
from various points of view. First we use a multipartite 
global entanglement as a measure. Secondly, we study 
the symmetric GHZ-like ground state and its transition 
to the paramagnetic product-state. Thirdly, by study- 
ing QPT in various dimensions we are able to establish 
common features of such a transition in different univer- 
sality classes. We find that global entanglement is a good 
indicator of such transitions with universal aspects, in- 
cluding scaling, in any dimension. The well-suitedness of 
such a measure is displayed in the nice fits obtained in 
figures such as Figure 2 or Figure 5, for example. As a 
by-product, we find critical points and various exponents 
for the Id and 2d models consistent with previous stud- 
ies. We note that our estimation of the critical points 
for the Id and 2d models are to within one percent of 
the generally accepted values, an impressive result given 
the limited size of the systems studied here, providing 
further evidence for well-suitedness of our measure when 
compared with similar studies using bi-partite measures. 
Our estimation of v, although acceptable, are under- 
standably less impressive as finite-size scaling collapses 
require larger system sizes to obtain better estimates for 
t/[47j]. We note that our main goal is to investigate the 
(universal) features of global entanglement in quantum 
phase transitions, and not to produce reliable exponents 
for such models. Since the parameter d determines the 
universality class of the systems considered here, the fact 
that we see similar behavior of global entanglement at 
the QPT regardless of d, shows what we have thus far 
referred to as universal features of global entanglement. 
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We close by mentioning that similar studies could be A. Langari for his comments on an earlier version of our 
carried out for more general spin models exhibiting more paper, 
complicated quantum phase transitions. It would be in- 
teresting to see if such universal features of global entan- 
glement carry over to other models. 
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